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ABSTRACT 

The  paper  addresses  the  problem  of  the  implementation  of  nonhomogeneous 
essential  Dirichlet  type  boundary  conditions  in  the  p-version  of  the  finite 


element  method. 
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1 .  INTRODUCTION 

^  A  recent  addition  to  the  field  of  finite  element  analysis  has  been  the 
development  of  the  p  and  Iv-p  versions  of  the  finite  element  method.  For  two- 
dimensional  problems  the  p  and  h-p  versions  have  been  implemented  in  the 
commercial  system  PROBE  by  Noetic  Tech.,  St.  Louis  with  a  first  release  in 
1985  and  the  second  in  1986'*£TJi'.  The  three-dimensional  commercial  finite 
element  code  FIESTA  having  some  p-version  capabilities  had  been  developed  at 
ISMES  (Institute  Sperimentali  Modelli  e  Strutture)  in  Bergamo,  Italy,  and  has 
been  available  since  1980  in  the  United  States.  An  implementation  of  the 
three-dimensional  p  and  h-p  versions  for  the  Cray  computers  is  presently  being 
done  at  the  Aeronautical  Research  Institute  of  Sweden  tFlygtekniska 

„  M  l.  ^  ■  .  /  'i 

Forsoksanstalten-FFA)  [1  ]  and  by  Noetic  Tech.  -  ■  }  > 

The  p  and  h-p  versions  are  being  used  very  successfully  today  in 
industry.  See,  for  example  [8].  A  survey  of  today’s  state  of  the  art  may  be 
found  in  [3]  where  pertinent  references  are  given.  For  basic  theoretical 
results  we  refer  to  [4],  [5],  [6],  [10],  [11].  For  some  implementational  and 
engineering  aspects,  we  refer,  for  example,  to  [8],  [14],  [15]. 

An  important  part  of  any  code  in  finite  element  analysis  is  the  problem 
of  imposing  nonhomogeneous  essential  boundary  conditions  of  Dirichlet  type. 

We  addressed  this  question,  among  others,  in  [5]  for  the  two  dimensional 
problem.  Our  technique  was  implemented  in  PROBE  -  Release  2,  and  tested  very 
successfully.  In  [5]  we  assumed  that  the  nonhomogeneous  boundary  conditions 
had  slightly  higher  smoothness  than  the  minimal  possible  one,  namely  we 
assumed  that  the  boundary  data  belongs  to  Hk(D,  k  >  1  (instead  of 
minimally  k  >  V2  ).  With  this  assumption  we  obtained  an  optimal  error 
estimate.  There  are  some  indications  that  for  %  <  k  <  1  the  method  proposed 


In  [5]  could  lead  bo  a  loss  in  the  convergence  rate  and  that  for 
k  <  3/**  there  exist  boundary  conditions  for  which  convergence  may  not  be 
achieved  at  all. 

This  paper  proposes  an  implementation  which  guarantees  convergence  for 
all  k  >  *4  .  We  get  the  almost  optimal  estimate 

|u  -  u  |  <  Cp"(k'4)  log%  |u|  y  . 

P  H1  (Q)  Hk  4(Q) 

In  addition,  we  discuss  some  other  possible  methods  for  imposing  the 
nonhomogeneous  boundary  conditions  of  Dirichlet  type.  In  a  forthcoming  paper 
we  will  discuss  general  boundary  conditions  with  partial  constraints  for 
systems  of  equations.  Such  conditions  are  important  in  the  theory  of 


elasticity,  for  example 
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2.1.  NOTATION 

Let  R2  denote  the  usual  Euclidean  space  with  x  -  (x1 ,x2)  €  R2.  By 

Q  c  R^  we  denote  a  bounded  Lipschitzian  domain  with  piecewise  smooth 

n 

boundary  38  •  r  *  U  T  where  r,  are  smooth  (open)  arcs.  The  end 

i  -  1 

points  A^,  i  *  1,...,n  of  iv  will  be  called  the  vertices  of  8.  For 
fixing  the  ideas  we  will  restrict  ourselves  to  simply  connected  Lipschitzian 
domains  although  our  results  hold  in  general.  We  will  also  consider  8  c 
R^ ,  i.e.,  a  »  I  -  (a,b) . 

Let  L2(8)  -  1^(8),  1^(8),  Hq( 8) ,  k  £  0  denote  the  usual  Sobolev 

spaces.  For  u  €  1^(8)  we  denote  by  |u|  .  and  lul  v  the  usual  norm 

Hk(8)  H 

and  seminorm,  respectively.  For  k  >  0  not  an  integer,  we  define  HK(8)  and 

|*|  .  by  the  K-method  of  the  theory  of  interpolation  ([9],  [12]) 

H*(8) 

H£*e(Q)  -  (H£(Q),  H£+1(8))0j2 

0  <  e  <  1,  t  «■  0  -  k,  1  integer. 


The  norm  is  defined  accordingly,  i.e.,  with 


(2.1a)  K(u,t)  -  inf  ( |v |  ,  ♦ 

v€H£(8),w€H1^  (8)  H 
v+w-u 


t  |w|  ) 

H*  '(8) 


we  define 

(2.1b)  M  -  (/  t~3K(u,t)2  & )\ 

IT  y( 8)  0  C 

Analogously  we  define  Hq(8)  -  ( Hq( 8) , Hq*1 ( 8) )Q  2  (k  a  noninteger)  with  the 
norm  |u|  . 

hJ  9(8) 
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We  will  be  especially  interested  in  the  one  dimensional  case. 


Let  Q  -  I  -  (-1,1).  Then  we  define  |u(  .  and  |u|  .  by 

fTd)  h£(i> 

(2.1b)  for  k  a  noninteger.  For  k  t  integer  +  ‘4,  we  have  for  u  e  Hq(I) 


(2.2) 


>  k 

H*(I) 


rti) 


where  by  »  we  denote  the  equivalency  of  norms. 


For  k  -  integer  +  %  we  have 


(2.3) 


duu,2 


M  M4v  -  (M  ♦  jo-*  >  o  ) 

H?  *  (I)  HW  ’  (I)  dxP  H  (I) 


This  space  is  often  denoted  by  Hqq  (I)  in  the  literature  and  we  will  use  this 


notation  as  well.  Analogous  results  are  true  for  general  Q. 


Mow  let  ft  be  the  space  of  2ir-periodie  functions.  For  u  €  H  we  will 


write 


(2.H) 


u(C)  -  l  a1  cos  j  5  ♦  I  b.  sin  j  £. 
j-0  J-1 


Then  we  define  Hk  and  |*|  by  (2.1a)  and  (2.1b)  as  before.  We  have  then 


(2.5) 


I  k  -  [  l  a2(1*j2)k  ♦  X  b2(1*J2)k]4. 

H  j-0  J  j-1  J 


Moreover,  if  a^,  bj  are  such  that  the  right  hand  side  of  (2.5)  is  finite, 
then  u  (  and  the  norm  |u |  is  defined  by  (2.5). 


For  I  an  open  Interval  or  straight  line  segment,  we  denote  by  s  the 


length  parameter.  Then  we  define  Hk(I),  Hq(I),  |u|  ,  |u|  ,  as 

Hi  I)  Hq(I) 

before  with  respect  to  s. 


So  far  we  considered  only  spaces  of  scalars.  We  define  now  (Hk(a))m 


with  u  €  0^(0))“,  u  •  (^,...,1^),  Uj  i  Hk(Q)  and 
m  2  % 

lul  v  m  “  (  l  lu«l  t,  )  and  analogously  |-|  ,  etc. 

( H  ( Q) )°  i-1  1  Hi Q)  (Hj(Q))m 

By  C^(Q),  respectively  C^(Q),  we  denote  the  space  of  all 

functions  with  continuous  derivatives  of  order  k  £  0  (k  integer)  on  Q, 

respectively  Q. 

The  set  of  all  algebraic  polynomials  of  degree  (total)  less  than  or 
equal  to  p  on  fl  will  be  denoted  by  Pp(fl).  By  Pp(Q)  we  will  denote  the 
set  of  all  polynomials  of  degree  less  than  or  equal  to  p  in  each  variable 

on  Q.  For  V  c  R2  a  straight  segment,  we  define  Pp(D  as  the  3et  of 

polynomials  on  r  of  degree  less  than  or  equal  to  p  in  s  and  by  Pp(r) 

we  denote  the  set  of  polynomials  vanishing  at  the  end  points  of  T. 

Let  I  -  (-1,1).  Then  we  will  deal  with  two  different  polynomial  bases 

on  I: 

a)  The  Chebyshev  polynomials  Tk(x)  -  cos(k  cos^-1^(x)),  k  - 

0,1 ,2 . 

b)  The  integrals  of  Legendre  polynomials  Pk(x), 

x 

^(x)  -  (2k-1)  \  Pk-1(t)dt  -  Pk(x)  -  Pk_2(x),  k  -  2,3,.. 

-1 

Obviously  the  3et  {TkJ,  k  -  0,1,.  ..,p  is  a  basis  of  Pp(I)  and  (^1,  k 
2 . p  is  a  basis  of  Pp(I). 


2.2.  THE  MODEL  PROBLEM 

Let  H  -  (h’U))™,  Hq  -  (  Hq  (  Q ) ) ra ,  and  let  B(u,v),  u  -  (u, . um)  € 

H,  v  -  (v^ . vm)  €  H  be  a  continuous  symmetric  bilinear  form  on  H  *  H 

satisfying 


»n»r 


r\cw 


V 
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(2.6)  B(u.u)  *  Y|u|2  ,  Y  >  0 

H*  (  Q) 

for  any  u  €  Hg.  Further,  let  F  be  a  continuous  linear  functional  on  H. 

Let  on  Ti  the  matrices  /^(s)  -  {a£*j[(s)},  k,fc  -  i  - 

1,...,n  of  smooth  functions  be  given.  We  shall  assume  that  the  matrices  A^ 
are  nonsingular.  Further,  let  g^  -  {g^^},  k  *  i  -  1,...,n  be 

a  function  defined  on  ri .  We  will  say  that  g  *  {g^l  €  (H^  (D)m  if  there 
is  a  U  €  ( H* ( Q)  )m  such  that  g  »  U|  p 

Let  now  h  -  {h^},  i  -  1,...,n  be  such  that  with  g^  » 

Aj1(s)h^i^,  g  ■  {g^1.^}  €  (H^(D)ra.  Then  our  problem  is: 

Find  Uq  <E  ( ( Q) )m  such  that 

а)  Voir.  *  hCl3  *  AigCl] 

(2.7)  1 

б)  B(u0,v)  -  F(v) ,  Vv  6  Hq  -  (H^B))™. 

Because  of  the  assumptions,  the  problem  has  a  unique  solution. 

Remark  2.1.  We  could  obviously  transform  a)  into  Uq|j,  -  A’^h^''  * 

1  i 

g^-1,  i  -  1,...,n  and  assume  that  Aj^  is  a  unit  matrix.  Nevertheless,  our 
formulation  is  more  general  and  computationally  natural.  For  example,  in 
elasticity  theory,  although  we  formulate  the  problem  in  displacement 
components  u,  v  (in  directions  of  the  axes  x^  ,  X2),  we  prescribe  the 
conditions  for  displacement  in  the  direction  of  the  normal  and  the  tangent  to 
the  boundary. 

From  the  general  theory  of  interpolated  spaces  and  our  assumptions 

Til  l/ 

about  Ai ,  we  obviously  have  hL  J  6  H  *  (r^). 

So  far  we  have  assumed  that  the  Dirichlet  boundary  conditions  are  given 
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on  the  entire  boundary.  Nevertheless  our  theory  can  be  easily  generalized 
when  on  one  part  of  the  boundary  un  is  not  constrained. 


2.3.  THE  p- VERS ION  OF  THE  FINITE  ELEMENT  METHOD 

Let  us  assume  that  domain  Q  has  been  partitioned  into  curved 

rectangles  and  triangles  (see  Fig.  2.1). 

9  _ 

Let  5  »  U  0.  where  8^  are  (open)  curved  quadrilaterals  or 
i-1 

triangles  called  elements  of  the  partition  of  Q.  The  vertices  of  are 

called  the  nodes  of  the  partition.  We  will  assume  that  the  vertices  of  8 
are  nodes  of  the  partition. 


Fig.  2.1.  The  scheme  of  the  partitioned  domain. 

By  S  -  (-1.1)2  and  T  -  U,n  |  0  <  n  <  U+1 )  /3,  -1  <  £  <  0;  0  <  n 
<  (1~0^3.  0  £  £  <  1}  (see  Fig.  2.2),  we  denote  the  standard  square  and 


rvA  A  A  A  A  A  -W 


I  onto  Tj  j,  where  the  relation  between  Fj^  and  Fifj  is  obvious  (realizing 
that  the  sides  of  T  and  S  have  length  2). 

Let  now  Pp(Q)  -  lu  €  ( Q)  |  ulF^1 (x^ , x2) )  €  Pp(T)  if  Qj  is  a 

curvilinear  triangle,  u(F ^ ; (x^ ,x2) )  €  Pp(S)  if  is  a  quadrilateral}. 

Further,  we  denote 


Pj(fi)  -  (Pp(C))m  and 
.-1 


( Up  Up>|n)  . 


Let  Up  €  P™(  £2) ,  then  up  k(F~  (x1  ,x2)) ,  k  -  1,2,...,m,  is  obviously  a 
polynomial  of  degree  p  on  every  side  of  T,  respectively  S.  We  will 
identify  the  sides  of  T,  respectively  S,  with  I  -  (-1,1)  in  the  obvious 

0  I/ 

way.  Assume  now  that  a  projection  mapping  Pp  is  given  which  maps  H  (I), 
k  >  l4  onto  Pp(I)  with  (PpU)(±1)  -  u(±1)  and  PpU  -  u  if  u  £  pp(  I).  We 
remark  that  Hk(I)  <^0(1)  for  k  >  \  and  hence  u(±1 )  is  well  defined. 

The  pr-version  of  the  finite  element  method  for  solving  our  model  problem 


consists  of  finding  up  £  ?”(«)  such  that 


a)  u  =  P  g  on  F. 

p  p  i 

(2.8) 

B)  B(u  ,v )  =  F(v ) ,  V  v  €  pm(Q)  n  (Hhfl))m. 

P  p  0 


Obviously  the  method  strongly  depends  on  the  choice  of  the  projection  Pp. 

This  choice  can  be  influenced  by  various  factors  like  accuracy,  implementa¬ 
tion,  the  type  of  problems  to  be  solved,  their  formulation,  etc.  For  example, 
the  matrices  in  the  model  formulation  in  Section  2.2  can  be  used  in 

(2.8)a  instead.  The  main  part  of  this  paper  consists  of  analyzing  some 


choices  of  this  projection. 
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.  THE  PROJECTIONS  Pp 
3.1.  THE  PROJECTION  Pj£  AND  ITS  BASIC  PROPERTIES 

o*4 

In  this  section  we  will  introduce  projection  operator  Pp ,  respectively 

l4»  A 

P  ’  ,  of  Hk(I),  respectively  (Hk(I))m,  onto  Pp(I),  respectively  pp(I). 

We  prove  first 

Lemma  3.1.  Let  I  »  (-1,1)  and  I  -  (O.tt)  -  F-1  C I)  where  F(S)  * 

4*  A 

cos  £,  £  €  I.  For  u(x),  x  €  I,  let  u(£),  (  (  I  be  such  that  u(cos  5)  ■■ 
u(0*  Then  for  any  u  €  Hk(I),  k  >  *4  ,  we  have 


(3-D 


(3.2) 


-  lui  l/  ^  C2I“I  y 

H4(  I)  I)  *  wh 


M  J.  >  C(»C)  |u|  ^  , 

H  (  I)  HK(I) 


where  0  <  C^  <  C2  and  C(k)  are  independent  of  u. 

Proof.  Let  Sq  -  [z  -  x+iy  |  jx|^1,  0  ^  y  <  1)  and  Q  »  U  -  £  +  in 

|  cos  C  €  Sq}.  (See  Fig.  3.1). 


z  plane 


z  =  cos  £ 


Fig.  3.1.  The  rectangle  Sn  and  its  image  Q. 
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Let  u  €  Then  there  exists  U  6  ^(Sq)  such  that 


M  1 


H  (SQ) 


<  C|u|  v  and  U| T  -  u.  Let  V(£)  -  U(cos  O-  Then 
H^I) 


and  also 


|  v|  1  -  |  u  |  , 

H '  ( Q )  H'(Sq) 


|V[  ,  <  C|U| 

H  (Q) 


H  <V 


Hence,  by  the  imbedding  theorem  we  have  with  v  »  V|-t 


M  V  .  -  M  u  «  <  C|V|  <  C|(j| 

H2(I)  H/2(I)  H*  (Q) 


£  cM  v 

H’(SQ)  H4(I) 


where  C  is  independent  of  u.  The  inequality 


M  v  ^  cl”l  Y  - 

H/2(I)  H^I) 


can  be  proven  analogously  by  changing  the  role  of  SQ  and  Q.  Hence  (3- 
proven. 

Let  us  prove  (3-2)  for  k  »  1.  We  have 


l^l2 

d*  h°(i) 


r  fdui2  7.  fdui2  /.  2 

l  (d5>  «  ’  I.  &  '  ’-x  d* 


-1 


«-1  .  2 
rdu> 


-<  I  ©d*“  @0 

-1  aX  QX*  H°(  I) 


-  2 
idui 


M  ,  .  <  1^1  Q  +  M  o 

H-  tTJ  ^  H°(  I)  H°(I) 


<  l^fo  .  *  c|G|  y  ^ 

04  H  (  I)  H/z(l) 


.~.2 


1 )  is 


Hence 


1  4 


<  M2t  +  C|u|2,  <  C|u|2 

H  (I)  H*(I)  H  (I) 


In  the  same  way  we  prove  the  inequality  for  k  >  1,  k  integer.  For 
k  2  l4  nonintegral  the  inequality  follows  from  the  standard  interpolation 
theorem. 


Let  u  €  Hk(I),  k  >  %  .  Then 


\  < 
—  V  1  2 


-1  (1-x&) 


and  we  can  write 
(3.3) 


u  ( x )  «  l  a.Tjx) 

j-0  J  J 


where  T^(x)  are  Chebyshev  polynomials. 


We  define 


i)  for  u  €  Hk(I),  k  *  \ 


(3.4a) 

PD2  U  ’  I 
p  j-0 

ajTj (x)  *  Pp(D 

ii)  for  u  € 

Hk(I),  k  >  \ 

(3.**t>) 

8*4  u  - 
P 

P*4  u  *  u 

where  u  is  linear 

on  I  and  (P^u)(±1)  -u(±l). 

Theorem  3*2. 

Let  u  €  Hk ( I ) . 

Then  for  k  1  V2 

(3.5a) 

Iu  ”  Pp  ul  ‘4 

<  C(k)p"(k'’4)|u| 

HK(I) 

and  for  k  >  \ 

p  H4(  I) 
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(3.5b)  |u  -  $  u|  £  C(k)p'“^k~l^|u|  . 

H^I)  Hk(I) 

(3.5c)  |u  -  u|  x.  £  C(k)p"(k“^log‘4p  |u| 

Hq2q(I)  HK(I) 

The  constant  C  in  (3-5a,b,c)  is  independent  of  u  and  p  but  depends  on 
k. 

Proof .  Let  0(0  -  u(cos  O.  Then  0  €  H  is  an  even  2»-periodic 
function  and 

m 

(3.6a)  0(0  -  l  .uf  cos 

1-0 


16 


and  hence 


|u  -  £*£u|  i  i  C(k)p"(k"^|u| 


fT(I) 


and  (3*5b)  is  proven. 


Let  us  prove  now  (3.5c)  for  k  >  *4  using  (2.3).  We  have  to  show  that 
♦1  \ 

(3.10)  [J  (u-Pj  u)2(1«-x2)-1dx]  i  Cp"(k^)log1^p  |u| 

-1  V  HK(I) 

which  is  equivalent  to  showing 


[J  (u-?1^  u)2(sin  0-1d5]  4  Cp-^k-^  log'^p  |u| 


Hk(I) 


Let  us  first  observe  that 


(3-12)  (u-$  u)  -  l  u2,(cos2j£-1) 

J-C*£] 


1  u21+1 (c°s(2j+1  )£  -  cos  t) 


where  by  [a]  we  denote  the  integral  part  of  a.  Hence 


»T  I  TT  « 

J  (u4£  u)2(sm  sr'os  S  2[f  (  I  u21(00S  2J5-D)  jjjpj  de 


0  j.[^3 


r  c  1 

J  (  l  U2J„  (cos(2J»\  )t-003  5)1  as] 


°  J-CEf] 


5  8[J  (  l  u2J  sln2jq  _J 


sin  £ 


I 

H  .  »  «  ■  «  •  .  *  .  • 
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tr  •  d  . 

*J(  l  u2J„  .inu*t)!  Sinn)  «i 


3in  £ 


,  .  e1(t)  e?u)  -  -1 
i  c/( (  I  |u2j|0.j2)^^ - SJlJlU) - J|ln 


2-e,  (0-e2(0 


0  **  (i+j2: 

k  e  (O  e_U)  2-e  U)-e?U) 

♦  (  X  K1+1|(i+j2)2  ^ 11  1 

j-tei]  J  - 


)  ]  rr~“F  ^ 


b+j2)2 


£  Cp'2(k'^  [J  - 

0 


2e1  (5)  2e_(0  2(e.  (£)+£.(€)) 

**  1  (IP-O  2  p  1  2 


sin  C 


d^iiur 


^(D 


sin  £ 


for  0  £  e,  (O  <  2~r!  -~  *  2e. 


Choose 


e.,(0  *  e  on  [0,  -],  •  0  otherwise 


e0(£)  -  e  on  [x  -  — ,  ir],  -0  otherwise. 

^  P 


Then  we  obtain 


1 

IT-  — 

P 


s  Cp'2(k'1^  [p2eJ  52e”!de  ♦  j  (Sin  ♦  p2eJ  (*~02e  1  dt ]|u |2 

0  1_  1  n  ( I) 

P 


IT-  — 

P 


S  Cp'2(k"^  log  p  |u| 


Hk(I) 


This  proves  (3-11  )  •  (3-5c)  follows  from  (3.5b)  and  (2.3). 


Remark  3  - 1  - 


The  necessity  of  the  term 


log%  in  (3.5c)  is  an  open 


question. 


Remark  3.2.  For  u  €  Hk(I)  with  k  >  %  ,  we  may  prove  (analogously  to 


(3.7)) 

(3-13) 


t—  —  «2  p  -2(k-t)  a  >2 

-  up|  <  Cp  |u| 

HX(I)  H  (I) 


%  <,  l  <  k 


Theorem  3.3»  Let  u£P  (I),  u(±1)  -  0.  Assume  that  |u|  i  A, 

P  H4(I) 


|u|  Q  <,  A.  Then 


CV(I) 


(3.1*0 


|u|  X,  i  CA  log*4  p. 

Hoo(I) 


Proof.  As  before  we  have 


C|] 


0  -  [  b21(cos2jt-1 )  *  l  b2i+1 (cos(2j+1  )e-cos  C) 

j-0  J  j-0  J  . 


and 


\  b?(1*J2)‘4  *  A. 

j-0 


To  prove  (3-14)  we  have  to  show  that 


it  ,2 

J  jrrr «  *  c  » *  ■ 


1 


We  have  analogously  as  before  for  e1  -  c  on  [0,  -  0  otherwise,  e 


e  on  [*  -  *],  -  0  otherwise,  e  >  0 


v  J2. 
— 
'  n 


ir~  — 


sin  £ 


-<  cllip«inE«M!o,?  *^£-4rlalV:  I 


£  C  log  p  A  . 


Let  us  generalize  the  operator  P^.  Assume  that  for  x  €  I, 


0  <  Sq  i  a(x)  £  <  •.  Let  u  €  Hk(I),  k  2  Then  we  define 

respectively  as  follows: 

for  k  i  %  : 


(3.15a) 


for  k  >  *4 


(3.15b) 


p,p(°p,p’°  u>  -  pp  ( au )  * 


p4,a  u  -  P^,a  u  +  u 


where  u  is  linear  and  (^,ot  u)(±1)  -  u(±1). 

We  will  show  that  (3-15a)  and  (3.15b)  uniquely  define  P^,a  and 


u(x)  -  l  u-.T.fx) 
j-0  J  J 


a(x)  -  l  a.T.U) 
j-0  J  J 


(3.16) 


U_  (  X  )  -  P^'a  u 

r  P 


1  Vj1 


Then  using  (3.6)  we  have 


(3.17a) 


I  u j  cos  K, 


(3.17b) 


5(0  -  l  a]  cos  j£. 

J-0  J 


(3.17c) 


ap(U 


l  bJ 


cos  JC 


where  {bj}  are  determined  from  the  system  of  p  ♦  1  linear  equations 

•  p 

(3.18)  \  [(I  o  cos  j£)(  l  b  cos  J5)]cos  ££  dC 

-i  j-0  J  j-0  J 


We  have 


’  H 

j  (  l  a.  cos  j0(  I  u  cos  J5)cos  U  d£,  l  *  0,1, 

-»  j-0  J  j-0  J 


(3-19) 


(  l  a  cos  j£)(  1  u.  cos  js)  -  I  c.  cos  j£ 
j-0  J  j-0  J  J-0  J 


where 


co  "  ‘4  C“ouo  ♦  l  Vk] 

k-0 


3J  -  \  [\  a  U  ♦  I  CL  U  *  l  <*k_jUk]. 

k-0  J  k-0  J  k-J  J 


By  (3.15),  (3.17c)  and  (3.18)  we  see  that  dp(C)  is  such  that 


(3.20) 


PP  auP  *  PP  au* 


Let  u  -  { <x  -  {a^}  ,  —  *  £  -  (ckJ  ,  k  *  0,1, ...  .Let  l 

m 

-  {a  -  {a,,}1,  k  -  0,1,...,  la  I,-  -  l  a|"  <  •} .  Since  5  is  bounded,  for 

~  2  0  K 

any  u.  €  t2  we  113 ve  £  €  l2  and  we  can  write 


(3-21) 


c  -  A(a)u 


where  A(a)  is  the  infinite  matrix  with  coefficients  stemming  from  (3-19). 
Matrix  A(o)  is  then  a  mapping  of  i2  into  t2.  For  any  integer  p  l  0  we 
denote  by  Ap(a)  the  (p+1)  *  (p*1 )  principal  3ubraatrix  of  A(a).  By 
(3.17c)  can  be  Identified  by  a  (p*1)  dimensional  vector  b  (  R*5*1  with 


(using  3-18) 


(3.22) 


Ap(a)b  -  Cc]p+1 


where 


{Cj}T. 


j  -  o. 


We  will  now  show  that  for  any  d  €  Rp+1 


<lT(Ap(a)d)  2  agCd^d) 


which  guarantees  unique  solvability  of  (3.22)  and  hence  the  existence  of 


pp2,a* 


Lemma  3-1*.  Let  u  *  {uj”  €  £->.  Then 

— —  —  J  C 

uT(A(a)u)  2  a0(uTu_)  -  aQ|u|^  . 

2 

m 

Proof.  Let  u  €  i2.  Then  0  -  l  u,  cos  €  L2(I),  (I  -  (-<■*,  it)) 

j-o 


J  (au)ud£  *  J  au^dS  >  a0|uj^ 


L_(  I) 


Using  this  with  (3-21),  we  get 


(3.23) 


T  If”  ’  U  ? 

utA(o)u  -  —  J  u(5Q)d£  2  —  |u| 

*  -¥  *  L2(I) 


which  yields  the  Lemma. 


Corollary  3.5.  Let  €  Rp+1.  Then 


aJ(Ap(a)ap)  -  aJutcOap)  2  50(ajap). 
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Corollary  3-5.  shows  that  Pp,a  u  exists  and  is  unique  when  Q  €  L^d)  and 
that 


(3.24) 


|aP  I 


pi  o  -  *  clal  o  . 

H°(  I)  HU(I) 


with  C  independent  of  0. 
We  now  note  that 


-  — 


I  (“jbi )sin  jC 
J-0  J 


(aa)(,>  - 

at 


l  (-jc, )sin  jt. 
j-0 


iCi) 


Defining  the  matrix  f '  '  and  the  vectors  u^\  c^  by 


<a(1))u  *  -J(A)U,  -  “Jbj .  (elV)4  -  '-jc,. 


’31 


,(D 


'if  ^  >i 


O) 


J  ’  -  '3 


we  3ee  that 


(3.25a) 


-J  (c_)  1  -  “J  ( Au) ,  -  ( A^ 1  }u) 


(3.25b) 


Q(1 )  -  •  I  (u'1  *).  sin  je 

j-0 


(3.25c) 
so  that 


(5a)(1)  “  l  (cu>)j3inj<i 

j-0 


(3.26) 


(u^))t(A^)u)  •  \  j  a(1  )  (5<3)(1  }d6 


-TT 


easily  verified  that 


M2 


Let  |au'(x)|  i  on  I.  Then  for  CQ  >  e  ♦  —  t  e  >  0,  it  may  be 

“50 
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j  a(1)(oa)(1)de  ♦  c0  J  a(5a)dc  >  c|aj2 

-1  ~w  h  ( X) 

with  C  >  0  independent  of  0.  Using  (3-23),  (3.26),  this  shows  that 

(u(1))T(A(1)u)  ♦  Cq(u)T(Au)  *  C{(u(1))Tu(1)  ♦  uTu}. 

For  u  «  (  RP*1,  this  gives 

(3.27)  <2p0  ,)T(A^1  Jap)  •  ColajHVp)  2  CKa*1  >  >Tap<1  >  •  aj  ^ ! . 

Using  the  relations  (3-21),  (3-22)  together  with  (3.25)  and  (3-27),  we 

obtain 


(3.28) 


Kl 


pV(i) 


<  C|u| 


H](i) 


Combining  (3*24),  (3.28)  we  get  the  following  theorem  by  a  standard 
interpolation  argument 

Theorem  3.6.  Let  0  <  aQ  £  a(x)  ^  a1  <  •>,  |o^  ^  (x)  |  £  M-, , 

u  €  H^I),  l  l  %.  Then  if  up  -  p'£’n  u,  for  0  S  '<  <,  min(1fi.). 


Kl  k  ,  < 

P  HK(I) 

C(k)  |u  |  _ 

«( I) 

• 

Theorem  3*7. 

Let 

u  €  Hk(  I)  ,  k 

>  ‘/j,  0  <  CtQ 

i  a(x)  i  a1  <  », 

|a(1)(x)|  *  M1  . 

Then 

for  k  >  l4 

(3.29a) 

[u  - 

pp'a  u«  y  < 

P  H/2(I) 

C(k)pf(k~,/2) 

M  k 
hk(i) 

and  for  k  >  V2 

(3.29b) 

|u  - 

01, 

pp  ul  y  S 

P  H4(I) 

C(k)p"(k_,/2) 

M  k 

HK(I) 

o 


■  *>»» 


i  t'l  4'«  i'a  |‘i  iw  it.  *,»  >a  >||  <|i  >.)  »,*  * j 


(with  u(±1 )  -  (P^'au)(±1)) 


(3- 29c) 


|u  -  M,a  u |  \,  <  C(k)p  (k  ^log^  M  k 

P  HqQ(  I)  “  Hk(I) 


Proof.  Let 


Since  P*4  u  €  Pp(I),  we  have 


and  we  see  that 


Up  -  p£’a  u 


01  -  u  -  P^  u 


n  4  ♦  Ct 

Olp  -  Pp»*  0) 


pl4.a(pV2  u)  •  p4  u 


«p  *  up  ~  Pp  u- 


By  Theorem  3.6  with  k  «  V2,  we  have 


l5pl  y  _  <  CI5I  v  - 

H/z(  I)  H  2(  I) 


which  gives 


\%  ~  Pd  ul  V  -  *  cl°  "  Pp2  u*  -4  - 

P  P  H/z(I)  P  H4(I) 


so  that  by  the  triangle  inequality  and  (3.1),  we  obtain 


|u  "  I  y  <  C|u  -  P2  u|  y 

P  H2(I)  P  H4(I) 

(3.5a)  now  yields  (3.29a). 

Now  let  k  >  *4.  Then  we  see  that  using  (3.8),  (3.15a),  for  V2  <  l  < 
min(k, 1 ) , 
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(3-30) 


| (au-oup)(C) |  *  j (a(u-up)  -  P^(a(u-up) ) ) (£) | 


*  Cp_(A'l/*)  |a(u-u  )  |  f 

p  H*(I) 


^  Cp-(l"^  |u-u  j  _ 

P  Hl(I) 


i  Cp"(^)  p“(k_£)  Ju| 

HK(I) 


using  (3-13).  From  this,  (3.29b)  follows.  (3-29c)  will  be  proven  analogously 
as  before.  We  have 


u  - 


P%'a  u 


l  d2i(cos  2j£  -  1)  +  l  d2i+i (cos(2j+1  )C~cos  O 
j-0  j  ”0 


wi  th 


(3-3la) 


l  d2(1+jV  <  Cp“2(k'U|ul2  , 
J  lT(I) 


k  >  H,  0  *  l  <,  1 


j-o 


and 

(3.31b) 


.  -  -(k-‘4) 


lu  -  P  u|  S  p  -  v  t-I 

M  H  (I) 


M 


1  k~V  1  i-1/ 

Let  St  -  min^  ♦  1 )  >  j,  z  -  .  Then  analogously  as  before  we  get 


1 


*  ‘  2 

J  p  (a-p‘4.“  .,'.2  — 1 

2 
P 


P  u/  sin  i 


,  _  ""(2k— 1).  ■  i2 

d£  <  Cp  .  log  p  |u| 

HK(  I) 


using  (3.31b).  Further 


1 


gl/„a  .2  1 


J  (a- P£’a  u) 
0 


P  ’  sin  £ 


d£ 


i  CP '2(R'l)  |  ^-'(j_ifL^)dc  |u|2 

0  0  (Id)  H  (  I) 
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i  Cp 


-2(k-«.)p-4e  |uj2  ,  Cp 
HK(I) 


-(2k-1)|u|2 


Hk(I) 


ii 

An  analogous  expression  holds  for  J  ^  •  d£.  This  proves  (3* 29c). 


Let  on  I  the  equation 


(3-32) 


au  -  f 


be  given  with  (£)  €  Hk(I),  k  >  l/2,  a(x)  >  aQ  >  0,  and  |a(1)(x)|  £ 

M1 .  Obviously  now  u  -  f/a  and  u  £  Hk(I).  Our  aim  is  to  find  up  £  Pp(  I) , 


(3.33) 


so  that 


(3.34a) 


respectively 


aP  “  2  a.Tjlx) 

P  j-0  J  J 


|u  “  u  |  <  Cp  U  4)|u| 


n  1  V 
P  H4(I) 


HK(I) 


(3.34b)  |u  •  u  |  v  <  Cp"{k"4)log^)  |u| 

P  H4  (I)  H1 


!k  (  I )  * 


We  have  seen  that  (3-34a,b)  can  be  achieved  so  that  up  »  P^(f/a), 
respectively  up  »  Pp,a  (f/a).  The  coefficients  aj  in  (3*33)  are  then 
determined  from  the  conditions 


(3.35a) 


respectively 


+1  +1  . 

j  -  u  T.  dx  -  j  — —  £  T.  dx,  j  -  0,1 . p 

p  J  *  J 


(3.35b) 


1  1 

f  -  au  T  .  d  x  »  f  -  f  T,  dx,  j  »  0 , 1  , . . . ,  p , 

-1  J~~2  P  J  -1  j  2  j 

/1-x  /1-x 


and  (3.34a)  is  achieved  for  k  2  %  .  For  k  >  %  we  achieve  (3. 34b)  so  that 


we  can  then  subtract  a  linear  function.  We  have  used  here  the  expansion 
(3*33).  Of  course,  we  can  interpret  (3.35)  as  a  projection  and  use  any  basis 
functions  which  are  proper  for  implementation.  For  example,  we  can  write 


u„ 


B 

l  biA 

k-0 


where  i|»0  «  1  ,  ^  -  x,  and  4>k,  k  >  2  are  as  defined  in  Section  2.1.  Then 

condition  (3.35a)  gets  the  form 


r  1  ,  .  f  1  f  .  . 

J  — —  up*k  dx  -  j  — : —  r  <*x, 


^  /-x2 


+1 

J 

-1 


a  k 


k  *  0v1,...,p 


'1 — x 


from  which  coefficients  bk  may  be  determined.  (3- 35b)  leads  to  an  analogous 
form.  Instead  of  (3.35a),  (3.35b)  we  can  also  use  the  transformed  form 


(3.36a) 

respectively 

(3.36b) 

where 


*1  «■»  j 

J  CL  cos  j?  dC  -  /  r  cos  jC  dt,  j  -  0,1 ... .  ,p 

-1  -ir  a 


*ir 


J  5G  cos  j£  dC  =  J  f  cos  jC  d£,  j  -  0,1,..., p 


-TT 


u  -  l  a.  cos  j£. 
P  j-0 


Fast  Fourier  transform  techniques  may  now  be  used  on  (3-36). 

The  mapping  P^’  °  can  be  generalized  easily.  Let  A(x)  -  {a,  .  (x)}  be 
a  positive  definite  n  *  n  matrix  (not’  necessarily  symmetric).  Let  u  » 
{u^,...^}  €  (Hk(I))ra,  k  2  l4*  Then  we  define  up  -  {up  m), 

UP.J  *  J  “  so  that 


28 


P-4(A  pj4,AUp)  -  P^(Au) 


where  P^  u  -  ( P^  u.j , . . . ,  P^  u,,,)  €  P®(I).  PpJ,A  is  defined  analogously. 


p  -  '-p  -i  p  -m'  *  p'*'*  'p 

Theorem  3.7  Immediately  generalizes  to 


Theorem  3.8.*  Let  |a[j^(x)|  <,  M-j .  Assume  further  that  yTAy  £  Qoly|2, 


aQ  >  0,  for  any  y  -  (y1,...,ym)  €  Rm,  |y |  -  £  y..  Then  for  k  £  % 

i-1 


(3.37a) 


|u-P^'Au|  y  . .  k 

P  (H/2(I))m  (Hk(I))m 


<  C(k)p"(k‘,/2)|u| 


and  for  k  >  \ 
(3.37b) 


|u-p£'Au|  y  i  C(k)p-(k-^  |u| 

(H/2(I))“  (^(1))” 


(with  u(±1 )  -  (P^*A  u)(±1  )) 


01. 


(3.37c)  |u  -  P^’A  u|  ,  m  <  C(k)p-lK“*'log*k>  |u| 
P  (Hl(I))m  ” 


-(k-V2),__l4 


‘00 


(Hk(I))m 


As  before,  we  are  interested  in  finding  up  so  that  it  is  an 


approximation  to  the  solution  of 


(3.38) 


Au 


f . 


Theorem  3.8  gives  a  constructive  way  for  the  determination  of  up  and 
provides  an  estimate  of  the  accuracy  obtained. 

We  have  assumed  in  Theorem  3.8  that  A  is  positive  definite  (although 
not  symmetric).  If  A  is  not  positive  definite  but  becomes  positive 
definite  after  permuting  columns,  than  such  an  approximation  having  the 

>/  A  i 

desired  properties  also  exists.  It  is  sufficient  to  use  Pp  instead 


of  pi*^  ,  where  A*  is  the  permuted,  positive  definite  matrix 


J 

.1 


.1 


The  positive  definitness  of  A  is  only  a  sufficient  condition  for 
existence  of  the  approximation  satisfying  (3. 37)  as  can  be  seen  from  the 
following  example.  Let  and  at22  be  con3tants  and  ^  i®22  “  a1 2a2 


>  OJq  >  0 

on  I.  Then  (3.38)  gets 

the  form 

(3.39a) 

o1 1  (x)u-, 

♦  a12(x)u2 

•m 

f1 

(3.39b) 

“21 U1 

*  a22u2  * 

f2 

and  hence 

(because  cu,.,  and  0^2 

are  constants) 

we  get 

(3-40) 

°21up  ,1  * 

a22up ,2  “ 

f2- 

Since  both  021  and  cannot  be  zero,  we  may  assume  that  *22  > 


(3.41) 


UP.2  ' 


'22 


and  hence 
(3.*»2) 


so  that 


((°11  - 


“l  2a21 


22 


)Un  ^ 

P.  < 


-  Pp  (*12PP  r2>  *  PP  ft 
22 


Using  the  fact 


ou  -  -i5-2S-  >  -!S-  >  -2-  >  0, 

°t22  °22  °22 


we  see  that  u  .  having  the  desired  properties  exists  by  Theorem  3.7. 
P  *  ' 


(3.41),  we  see  that  u_  o  also  has  the  desired  properties. 
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3-2.  NUMERICAL  ASPECTS  OF  P^ 

Let  us  consider  the  problem  (3-38)  with  m  -  2  on  I  -  (-1,1) 


a1l 

U1  *  a1 2U2 

•  f1 

a21 

U1  *  a22u2 

»  f2 

Assume  that  f-, 

and  f2  are  selected  so  that  u 

■)  »  |  x  | 1  * 7  and  u,  » 

cosh  x  | x | 1  * 7 . 

We  will  now 

consider  coefficients 

a,  ,  of  various  smoothness 
*  >  J 

I) 

- 

sin  x,  a12 

• 

cos  x 

a21 

■ 

-cos  x,  au,2 

M 

sin  x 

II) 

4?J 

m 

x^2 

for 

x  <  0 

xq+1 

for 

x  2  0 

4S1 

m 

(X  ♦ 

for 

x  < 

(x  ♦  %)*" 

for 

x  >  -l4 

r  a  i 

n  uM  J 
a21 

m 

(X  -  '4)0 H 

for 

x  <  l4 

(x  -  V*)^2 

for 

x  >  l4 

m 

-x*" 

for 

x  <  0 

-xq+2 

for 

x  >  0 

III) 

a11 

- 

1 

for 

-i  <  X  ^  -‘4 

0 

for 

-l4  <  x  ^  o 

1 

for 

0  <  x  <,  % 

0 

for 

V,  <  x  <  1 

0 


for  -1  <  x  £  -*4 


*12 


1 

for 

-l4  <  x  i  o 

0 

for 

0  <  x  i  '( 

1 

for 

X  <  X  <  1 

a21 

-  0 

for 

-1  <  X  <  -  3/“ 

2 

for 

-  3/H  <  x  i  0 

0 

for 

0  <  x  1  ‘4 

2 

for 

X  <  x  <  1 

a22 

-  -1 

for 

-i  <  X  £  -‘4 

0 

for 

-X  <  x  1  0 

for 

0  <  X  4  ‘4 

0 

for 

l4  <  x  <  i . 

-  PV> 
P 

,  and  n.(p)  = 

2f  T\  1 

lu.  -  U.  I  1, 

1  1  I’P'uV 

We  compute  c^ip)  *  j 

“  v  y  2  “  ' 

where  ut  is  the  ith  component  of  P'*’  u.  Obviously,  £  nr  The 
purpose  of  the  computation  is  to  see  the  effect  of  the  matrix  A  on  the  error 
of  the  approximation,  especially  in  dependence  on  the  smoothness  of  the 
coefficients  a^.  Table  3.1  shows  some  of  the  results.  We  see  that  the 
influence  of  the  smoothness  on  the  performance  is  small,  as  expected. 

Moreover,  we  see  that  the  asymptotic  rate  of  convergence  from  Theorem  3.8  is 


$ 

> 


i 


Vl 


already  observed  for  small  p. 


32 


TABLE  3.1 


The  errors  e1  and 


1 


t 


. 


* 

i 

! 


p 

C1 

n1 

e2 

o2 

CASE  I 

4 

2.2683-2 

2.2718-2 

2.1070-2 

2.1150-2 

8 

8.2301-3 

8.2586-3 

7.8898-3 

7.9168-3 

16 

2.7841-3 

2.791 0-3 

2.7534-3 

2.7601-3 

32 

9.1103-4 

9.1 242-4 

9.0843-4 

9.0977*4 

CASE  II,  q 

-  3 

4 

2.2683-2 

2.3788-2 

2.1070-2 

2.4886-2 

8 

8.2301-3 

8.6852-3 

7.8898-3 

8.9047-3 

16 

2.7841-3 

2.9217-3 

2.7534-3 

3-1241 -3 

32 

9.1098-4 

9.4950-4 

9-0843-4 

9.9006-4 

CASE  II, 

q  ■  0 

4 

2.2683-2 

2.3H7-2 

2.1070-2 

2.363 2-2 

8 

8.2301-3 

8.3735-3 

7.8898-3 

3.7462-3 

16 

2.7341-3 

2.8326-3 

2.7534-3 

2.9975-3 

32 

9.1098-4 

9.2162-4 

9.0849-4 

9.5976-4 

CASE 

III 

4 

2.2683-2 

2.2763-2 

2.1070-2 

2.1459-2 

8 

8.2302-3 

8.2438-3 

7.8898-3 

8.1540-3 

16 

2.7841-3 

2.8078-2 

2.7534-3 

2.9632-3 

32 

9.1 104-4 

9.2687-4 

9.0843-4 

1 .0191-3 

3.3.  the  PROJECTIONS  £p  AND  P°  AND  THEIR  PROPERTIES 

•*/ 

In  the  previous  section  we  have  been  interested  in  the  projection  Pp2 
We  can  introduce  the  ?p  projection  as  follows.  Let  u  €  Hq(I),  then 

P 

( 3 « **3)  Plu  -  I  a  *  (x) 

v  k-2 

where  iiik(x)  were  defined  in  Section  2  and  the  coefficients  ak  in  f  3  - 14  3 ) 
are  determined  from  the  conditions 

♦  1  *1 

( 3 .  uu )  f  (f>'  u)’4»]  d  x  •  J  u>'  dx 

-1  P  J  -1  J 

♦  l 

-  -  J  uK»!J  dx,  J  -  2,  .  . 

-1  J 

(where  the  3pan  of  {*j}  is  the  set  of  all  polynomials  or  degree  p  -  2). 
(3.1*1*)  may  be  used  to  define  Ppu  for  u  €  H°(l)  as  well. 

If  u  €  H1 ( I) ,  then  we  define  analogously  as  before 

( 3-**5)  £pu  -  v  «•  ^(u-v) 

where 

v  -  j(-x*1)u(-1)  ♦  (x*1  )uO  ) 

and  hence  Pp  is  well  defined  for  all  u  €  Hk(I),  k  >  V2.  In  [5]  we  have 
proven 


Theorem  3-9.  Let  u  (  HK(I),  k  >  1.  Then 


(3. U6) 


lu  -Ppul  v  <  Cp  U  4)|u| 

Hj0(I)  H  (I) 


We  see  that  this  theorem  13  analogous  to  Theorem  3-2  (in  fact,  it  gives  a 
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slightly  better  estimate).  Nevertheless,  Theorem  3*9  assumes  that  k  >  1.  We 
conjecture  that  for  %  <  k  <  1 


(3-47) 


3UPk  H  k 

u€Hq(I)  3(1) 


*  Pjuj  >  Cp 

W  H*(I) 


(Izi) 

(k-VJ.  2  ' 


In  [2],  [7]  we  have  proven 


Theorem  3 . 1  0 


s;,p  -pT .  lu  •  ^“1 

u€Hq(I)  HJ(I) 


P  1  0 
H  (I) 


i/ 

2  Cp*. 


|u-  fju|  S  Cp'4  |U| 

H  !I)  H  i!) 


From  Theorem  3.10  we  get 


(3-48) 


y  <  Cp4  |u|  , 

«00(I)  H00(I) 


but  not  necessarily  (3.47).  Nevertheless,  numerical  experimentation  suggests 
that  (3-47)  colds  for  k  » 

To  snow  it  let  for  a  >  0 


(3.49) 


u3(x) 


100 

[  dj^tx) 

i-1 


Select  a  -  -  and  compute 
P 
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XP<«> 


|u  -  p’u  I  , 
■  a  p  a*  V 


H^I) 


np(a) 


TC  xP(o)- 


Table  3-2  shows  the  values  of  n_  (a),  a  »  —  which  clearly  indicate  our 

P  p 


conjecture 


TABLE  3.2 


The 

values  of 

np(a), 

2 

a  -  -  . 

P 

P 

nP 

P 

np 

1 

•  371 

1 1 

.645 

2 

.479 

1  2 

.647 

3 

.539 

13 

.647 

4 

.576 

1  4 

.648 

5 

.600 

15 

.648 

6 

.616 

16 

.648 

7 

.626 

17 

.668 

8 

.  63^ 

18 

.648 

9 

.639 

19 

.648 

10 

.643 

20 

.647 

So  far  we  considered  the  projection  Pp .  Let  us  briefly  also  consider 
°0 

the  projection  Pp>  which  is  given  once  more  by  (3.43)  but  with  ( 3 .  ti U ) 
replaced  by 


(3.50) 


T  ~  I 

j  (  PpU )  4ij  (x  )dx  -  J  u4ij(x)dx. 


Using  the  results  from  [2],  [7]  it  is  possible  to  prove 
Theorem  3  •  11 

sup  FT  ,  Iu  "  11 1  i  2 


u€  Hq(  I)  H0(I) 


^p  11 1  •)  i  CP‘4 
h'(D 


W»vWV.**\* 


ViM1  y^iu^  nviuiwi  v^\rw.  iTKi/wwiuiu » 
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|u  -  pju|  <  Cp^|u| 

r  ,.'  /  t  \  «»' 


‘4. 


H0(I)  H0{1) 


This  theorem  suggest  that  (3.^7)  holds  also  for  the  projection  P;: 


3. 4.  COMPARISON  OF  THE  PROJECTIONS  Pp,  Pp  and  Pp 


We  are  mainly  interested  in  the  approximation  properties  with  respect  to 


the  norm  |*|  x.  .  The  previous  results  show  that  seemingly  the  most  robust 
H00(I)  v 

o  *2 

projection  is  the  projection  Pp  which  leads  to  the  error 


(3-51 ) 


2*4 


|u  -  P  u|  y  <  C(k)p  (k  /z)log^  |u| 

“  ti'2  /  r  \  n* 


"oo'1’ 


H  ( I) 


for  all  ‘4  <  k.  (We  conjecture  that  the  term  log4  p  is  not  needed.) 


The  projection  Pp  leads  to  (3-51)  also  (in  fact,  without  the  term 


V  13 

log7*  p)  only  for  k  >  i.  For  —  <  k  <  j-  the  possibility  exists  that  the 


projection  Pp  may  not  converge  at  a)  1.  In  addition,  the  projection  Pp*  can 


easily  be  generalized  to  p4*A  .  This  indicates  that  the  projection  p4  is 


preferable,  nevertheless  the  projection  Pp  is  almost  as  good.  On  the  other 


hand,  in  context  of  the  implementation  in  a  finite  element  code  using  shape 
functions  based  on  i(jj(x)  as  in  the  code  PROBE,  the  projection  Pp  is 


slightly  preferable.  The  projection  Pp  seems  to  have  no  advantages. 


4.  THE  p- VERSION  OF  THE  FINITE  ELEMENT  METHOD 

Let  us  consider  here  the  convergence  of  the  p-version  of  the  finite 
element  method  when  Dirichlet  boundary  conditions  are  prescribed  as  in  (2.7), 
(2.8).  We  shall  prove  the  following  theorem 


Theorem  4. 1  . 

O 


the  projection  Pp 
(4.1)  |u  -  g[l] 


Let  the  p-version  of  the  finite  element  method  be  based  on 
(see  (2.8))  such  that 


m  - 


<  C(k)f(p,i)p 


•(k- 


"Vl]i 


k-1/  r 

or  2(r.  ))• 

i 


and  let  the  exact  solution 


(4.2)  |u0 


u 

P 


( H1 ( Q) )m 


uQ  €  (Hk(«))m,  k  >  1.  Then 

C(k)p  1  ^maxCf(p.i) ,  log/z  p]p  ^  1  ^  Ju  |  , 

i  0  (HK(3))m 


Proof.  For  simplicity  of  notation,  we  let  m  ■  1  .  Let  Q  -  T  or  S. 

Let  1)^(5, n)  be  defined  on  Q  so  that  1^(5, n)  -  Ug(Fi(L,n)).  Because 

is  assumed  to  be  smooth,  U,  €  Hk(Q).  Hence  as  in  [5]  there  exists  U,  € 

1  1  t  P 

Pp(T),  respectively  Pp(S),  such  that 

(4.3a)  U , ( A • )  =  U,  ( Ai ) 

A  J  1  *  K  J 


where  {A.}  are  the  vertices  of  Q 


On  every  rifj  t  T.  Ut >p(Fi  (x ,y ) )  -  Uj fP(Fj  (x ,y ) )  -  ^ >p  -  Uj fP  -  <Pt , j ,p  * 
and  (see  Section  2.3  for  notation)  *i,j,p^  "^i.j.P  5  (  I  is  a 

polynomial  of  degree  p  in  one  variable  with  4^  j  p ( ±1 )  *  O.By  the  imbedding 
theorem 


(4.4) 


l*i  j  pi  v  <  cP“(k"1)[|u  |  -  |u  | 

i,J,P  H2(  I)  H  (Q)  J  H*(Q) 


<  Cp'^“1}  |u  | 

H  (Q) 


Using  Theorem  3-3.  (4.3c)  and  (4.4)  yield 


(4.5) 


l^i  j  pi  y  <  Cp'U'1)log4p  |u  1  . 

1,J,PHq20(I)  HK(Q) 


Applying  now  Lemma  4.7  of  [6]  there  is  a  vi  j  e  Pp(T),  respectively 


Pz  (S) ,  such  that 
P 


ri.j  "  vij(FI1(x*y))  *  n.j.p  on  ri.j* 

'i,J  *  0  on  3ni  •  ri,j 


lvi  j  I  1  *  I^i  j  pi  y  <  Cp“(k’°logS  |u  I  . 

l,J  H(n.)  1,J,P  Hq20(i)  u  hk( n) 

Let  now  r.  c  3Q  -  r  and  let  r5  c  .  Denote  gj;^  €  Pn(r{)  such  thai 

~  r„ 

gp  -  ur,p  “  Uq  -  up  +  gp  -  g 


and  hence 


lgD  “  ur  DI  V  -  Cp_(k  1)[1°8/2  P  +  f(p.t)]IuQ| 

P  ,P  h1(tj  U 


Hk(  fl) 


0  as 


and  <PrfpU)  -  (gp-ur<p)(Fr(0)  satisfies  tpn>p  €  Pp(I).  4>r,p(±1)  ‘ 
before.  Hence,  in  the  same  way  as  before,  by  Lemma  4.7  of  [6]  we  can 
construct  vr>i  such  that  vr  £  “  sp  ~  ur  p  on  r£*  vr  £  "  ®  on  ^ar 

r£  and  lvr  j, I  i  <  ISD  "  ur  J  u 

Pf  h'(Q_)  P  Ho0(r£) 


Defining  up  €  Pp(fl)  fl  H  ( fl)  by 


(*1.6) 


ui  “  l  vi.j> 
j 


we  get 


( *4.7) 


lun  “  QdI  1  <  C(k)p  (k  ^maxClog^p  +  f(p,i)]Ju-| 

PH  (fl)  i  °  Hk(n) 


Op  -  gp  on  r. 


This  leads  immediately  to  (4.2).  In  fact,  by  (2.8) 


B(up,v)  -  B(u0, v) , 


and  hence 


B(VVv)  "  B(u0  '  °p’v)’ 


V v  €  P p(Q)  0  hJ(C) 


Vv  6  p  iq)  n  hJ(o) 


where  up  is  defined  by  (4.6).  Because  up  -  up  €  Pp(fl)  fl  Hq(A),  let  v  * 


up  -  up.  Using  (2.6),  we  get  immediately 


luP  -  aplHi(n)  *  cK"apii 


H  (fl) 


and  hence  also 


|u0  ‘  “pi  1...  -  |u 


H  ( fl) 


‘JD  I  1  +  lUQ  ~  l"*p  I  •) 

P  H  (fl)  P  H  (fl) 


<  C|uQ  -  a  | 


H  (fl) 


WKWiWWW 


no 


This  gives  (4.2).  a 

Theorem  4.1  leads  to  the  following  corollary 

Corollary  4.2.  Let  Uq  €  Hk(Q),  k  >  1,  be  the  solution  of  the  problem 

(2.7),  Uq |  p  -  i  =  1,...,n.  Let  the  p-version  of  the  finite  element 

i  V  V 

O  <2  o  '2  »  i 

be  used  with  Pp  -  Pp  or  Pp  =  Pp  where  are  positive  definite 

matrices  with  C1  coefficients  on  (see  (2.7),  (2.8)).  Then 

lup  '  uol,..t 


<  C  P~<l‘"’  ’log'^  |u  I  . 

(H’(Q))m  (Hk(n))m 


Remark  1 .  A  similar  estimate  follows  from  Theorem  4.1  for  the  projec¬ 
tion  Pp.  In  fact,  for  this  case,  the  following  result,  proved  in  [53,  is 
true. 

Let  Uq  £  Hk(fl),  k  >  3/2  be  the  solution  of  the  problem  (2.7),  uQ | r 
*  g^*^,  i  *  1,...,n.  Let  the  p-version  of  the  finite  element  be  used  with 
Pp  -  °pl-  Then 


*V ‘ 


,  _  -(k-1 ) ,  , 

<  Cp  |uQ| 


( Hk(  n)  )r 
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